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Abstract 

O ■ Thermodynamic stability of composite vortex in a two-component superconductor is investigated by the 

Ginzburg-Landau theory. The predicted nature of these vortices has recently attracted much attention. Here 

we consider axially symmetric quantized vortex and show that the stability of vortex depends on three 

CN ■ independent dimensionless parameters: ki, K2, k^, where /c,(/ - 1,2) is the Ginzburg-Landau parameter of 

r^ , individual component, k^ = ^1/^2 is the ratio of two coherence lengths. We also show that there exists 

o 

O i thermodynamically stable vortex in type-l-i-type-2 or type-2-i-type-2 materials over a range of these three 

in 

0_i' parameters. 
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PACS numbers: 74.25.Ha, 74.25.0p, 74.20.De 



The existence of quantized vortex in a type-2 superconductor (SC) is one of the most striking 
phenomena in condensed matter physics [1]. The criterion for stability of vortex in a conventional 
SC is the Ginzburg-Landau parameter k, which is defined as the ratio of penetration depth to coher- 
ence length: k = A/^ [2]. Vortex can exist as a thermodynamically favorable state under external 
field in a type-2 material with k > 1/ V2, while the penetration of vortex is not thermodynamically 
favorable in a type-1 SC with k < 1/ V2. 

Recently, there has been growing interest in investigating the vortex in multi-component SCs 
ISl-lsl]. The predicted nature of vortex in these materials is quite different from that of vortex 
in a conventional type-2 SC. Babaev and Speight showed that interaction potential between two 
vortices in a two-component SC can be non-mono tonic: intervortex force is attractive at long 
range and repulsive at short range 151]. The key question regarding Babaev and Speight's work is 
whether vortices in a two-component SC are thermodynamically stable or not. In this paper, we 
revisit this issue addressed in Ref. [5]. The idea is that, vortex can survive as a thermodynamically 
favorable state if the Gibbs energy of the vortex state under the thermodynamic critical field is 
smaller than that of the fully superconducting state (Meissner state). We find that, the stability 
of vortex depends on three independent dimensionless parameters: k\,K2,k^, where Ki{i = 1, 2) is 
the Ginzburg-Landau parameter of individual component, k^ = ^1/^2 is the ratio of two coherence 
lengths. We also find that vortex is thermodynamically stable in a type-l-l-type-2 or type-2 -I- type-2 
SC over a range of these three parameters. 

Based on the Ginzburg-Landau model, we use free energy density in the two-component SC as 
follows: 

/ = U + Yj ^K^ - ^^^'^'l' + ^(1*^1.21') + nCvyi + ^i^P + ^(V X A)\ (1) 

;=1 ' 

where /„o is the free energy density of the body in the normal state in the absence of the magnetic 
field, Vd*?,!^) = a,|*F,P -I- b,-\^if/2 (i = 1,2). 77 is a coefficient characterizes Josephson cou- 
pling between two superconducting components. In the following we consider in particular weak 
Josephson coupling limit and set 77 = 0. We also assume that the effective mass m* and charge e* 
of two components are equal: m* = m*, e* = e*. There are four characteristic lengths: the penetra- 
tion depth Ai and coherence length ^; for each component are given by: At = {m* c^ / Ane*^'^^^)^^^ , 
^j = %/{2m*\ai\y^^, where ¥,0 = {-Uilbi)^^^. The thermodynamic critical magnetic field of the 



individual component is Hct(i) = Oo/(2 y/lnAi^j), where Oq = hc/e* is the flux quantum. The mag- 
netic field penetration depth and the thermodynamic critical magnetic field of the system ([T]) are: 
A = a/A] + \IAlr'l\ //,, = (//2^i) + Hl^^^yi\ Note that A < min(ii, ^2), H,, > max(//,,(i), //..p))- 

We consider axially symmetric quantized vortex in the model Q: 

*Pi = I^Fil e'\ ^'2 = |*P2l e'\ A = A(r)ee. (2) 

In order to study the stability of this vortex, we consider the Gibbs energy difl'erence between 
the vortex state under the thermodynamic critical field Hd and the fully superconducting state 
(Meissner state, or fully normal state under the thermodynamic critical field since these must be 
equal) 

A G = Gyortex(Hct) - Gm = Gvortex{Hct) - GniHct). (3) 

Let us show that if AG < then the isolated vortex can appear as a thermodynamically favorable 
state under external field H < Hct- There are three possible states: fully superconducting state 
(Meissner state), vortex state and fully normal state. And we note both the Gibbs energies of vortex 
state and normal state are decreasing with increasing field. If AG > 0, the normal state becomes 
energetically favorable when the field exceeds the thermodynamic critical value and vortex state 
is energetically unfavorable. If AG < 0, vortex state becomes favorable under certain value of the 
external field which is smaller than the thermodynamic critical field H < Hct. 

The Gibbs energies of the vortex state and the Meissner state can be written as 

GvortexiHct) = I dvgy, Gm = Gn{Hct) = I drgM. (4) 

Here ^v = /,, - Her • ( V x A) /4n,gM = /m = fno - H], / %n are the Gibbs energy densities of the vortex 

state and the Meissner state, respectively. To investigate the vortex stability in the considered case, 
one can write down 



A 



G = J Jr I ^ ^ |(V - ^A)¥,|2 + y(|¥i,2h + ^(H,, - V X A)4 . (5) 



We shall use instead of the variable r, the functions ^/ and A the dimensionless quantities 
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In the following we calculate AG and find 







(6) 



(7) 

Here Aj = 1//^^ + ^2/^2^ g^ ^ (^2 + K^Kp/iKl + K^Kp, Ci = KJ/iKJ + K^Kp, A2 = l/KJ + 

IIk\k], B2 = (KJ + k\kPI[k]{kI + K\Kp], C2 = K^KpiKl + K^Kp, Ki = Aij^i = 1,2) is the 
Ginzburg-Landau parameter of individual component, k^ = ^i /^2 is the ratio of two coherence 
lengths. The Ginzburg-Landau equations which determine the profile of the vortex solution are 
determined by minimizing the AG with respect to functions (A,(z =1,2) and A 

Alii, + — -%+ y[2MBi—A - \b,A^iIj, + i/r, - ,^1 = 0, 
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A2(«A2 + — - ^) + ^2A2B2—A - ^52AV2 + «A2 - <A2 = 0, 

P P- P 2 



A' A _ ,2 ^ ,2.. .^ /2AiiA? ^ /2A7iA2 



A' + - - ^ = (Ci<A? + C2<A^)A -(CiJ—^^ + C2 J-^-y (8) 

p p^ y Bi p y B2 p 

The order parameters vanish in the vortex core which is the phase singularity of both components. 
And it can be expected that the potential A correlates linearly with p in the core region due to the 
constant value of the field inside the core. Then the boundary conditions at p = are 



^i(0) = <A2(0)=A(0) = 0. (9) 

We assume the following power series solutions to Eqs. ([8]) in the region of the core: 1/^1 (p) = 
i;r=o bnP", ^lip) = Zr=o Cnp", A = i;r=o ^np" and can prove that 

'Ai = ^iP - vr( ^I^AiByai + \)bip' + Oip'), (A2 = cip - — -( yl2A^_ai + l)cip^ + 0(p^), 



A = aip - i(Ci ^J^bl + C2 ^J^cpp' + 0(p'). (10) 

Here a\,bi,Ci are three constants which will be deduced from the solutions far from the vortex 
core [lOJ: 
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FIG. 1: Stable vortex solutions in a two-component SC. (a) is a vortex solution in a type-l-i-type-2 SC, while 
(b) and (c) are vortex solutions in the type-2-i-type-2 SCs. Note the similarity of the vortex configurations 
in (a) and (b). 

TABLE I: Stable vortex solutions in a two-component SC 



Kl 


K2 


Kf 


bi 


C[ 


ai 


AG/iH^,A^/4) 


3 


0.5 


1/12 


1.7253 


0.4111 


0.2830 


-0.9044 


3 


3 


1/6 


1.8973 


0.5843 


0.2762 


-0.5721 


3 


3 


0.7 


1.6275 


1.2383 


0.3425 


-1.1159 



1 /I 

^i(^) = M^) = 1, B(cx,) = (-_(pA))U = 0, 

pdp 



(11) 



where B = (l/p)d(Ap)/dp is the scaled magnetic field: 5 = |V x A| ///,.,. It is clear from ©, ([8]), 
dH) and (ITT)) that the sign of AG, and the stability of vortex in the model are determined by the 
three independent dimensionless parameters: ki,K2,k^. 

We have numerically solved Ginzburg-Landau eqs. dH) with boundary conditions (fTOl) . (fTT]) 
and identified thermodynamically stable vortex solutions in the type-l-i-type-2 and type-2-l-type-2 
SCs, as predicted previously by the surface energy calculations [llj. Figure 1 illustrates several 
examples of stable vortex solutions. We found that the vortex stability in a two-component SC 
depends not only on the Ginzburg-Landau parameter of individual component, but also on the 
third parameter k^ = ^1/^2- 

Here we want to relate works in Ref. p] to the results of the present work. In their paper, 
Babaev and Speight identified three characteristic lengths in the model ([U): penetration depth A 
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FIG. 2: Examples of stable vortex solutions with ^i/A - V^/S and ^2/^ = 4 V2. (a) is a vortex solution in 
a type- 1 +type-2 SC, while (b) and (c) are vortex solutions in the type-2+type-2 SCs. Detailed parameters 
can be found in the table 2 below. 



TABLE II: Stable vortex solutions with ^JA ^ V2/3 and ^2//* ^ 4 V2 



Kl 


K2 


K^ 


fiM 


6/^ 


^1 


Cl 


fli 


AG/{HlA^/4) 


2.3 


0.4574 


1/12 


V2/3 


4V2 


1.4096 


0.3753 


0.3161 


-0.9122 


2.15 


1.0859 


1/12 


V2/3 


4V2 


1.4194 


0.3813 


0.3197 


-0.6696 


2.13 


1.9602 


1/12 


V2/3 


4V2 


1.4208 


0.3821 


0.3203 


-0.6387 



and coherence lengths of two components ^1,^2- They presented an example of vortex solution 
with ^i/A = V2/3, ^2/^ = 4 V2. Our results reveal that the intrinsic parameters which determine 
the magnetic properties of a two-component SC are three independent dimensionless parameters: 
Ki,K2,K^. It is easy to verify that the values of these three parameters are not unique for given ^JA 
and ^2 //I- This means that, the vortex solution presented in Ref. [5] is not the unique solution. In 
figure 2 we present several examples in which k\ ,K2 are different while ^1 /A = V2/3, ^2/-^ = 4 V2. 
It is clear from table 2 that the slopes of i/r 1,1^2 and A near the center of vortex core are slightly 
different. This proved that these vortex solutions are different solutions of Ginzburg-Landau eqs. 
([8]). In particular, there are cases in which two component are both of type-2 ((b) and (c) in the 
figure 2), which were not mentioned in Ref. [|5D. 

A notable aspect of our results is that in a type- 1 -I- type- 2 SC, stable vortex solution always 
has a extended core associated with the type-1 component (see fig. 1(a), fig. 2(a)). It is pointed that 
non-monotonic interaction between vortices originates from this exceptional vortex configuration 



||5i lZD. However, in a type-2+type-2 SC, the situation is much more complicated. We found 
that, stable vortex solutions in a type-2+type-2 SC may have a extended core (whose range is 
much larger than that of penetration depth of the system, see fig. 1 (b),fig.2(b),2(c)), or may have 
contracting core (whose range is smaller than the penetration depth, see fig. 1(c)). This means that, 
in a type-2+type-2 material, intervortex fore may be attractive at long range and repulsive at short 
range, as that of vortices in a type- 1 +type-2 SC. Alternatively, intervortex force in a type-2+type-2 
SC may be repulsive at all range, as that of vortices in a conventional type-2 SC. 

In conclusion, we have identified the intrinsic parameters which determine the stability of vor- 
tex in a two-component SC. Isolated vortex can appear as a thermodynamically stable state in a 
type-l-l-type-2 or type-2 -l- type-2 SC over a range of these three parameters: Ky,K2, k^. 
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